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1 .  INTRODUCTION 

At  present  the  designer  of  small  aircraft  must  insure  that 
the  structure  be  capable  of  withstanding  the  load  produced  by  a 
single  gust  of  specified  form.  However,  the  continuous  and  random 
nature  of  atmospheric  turbulence  suggests  that  the  designer  con- 
sider working  with  histograms  of  the  major  loads,  i.e.,  the  number 
of  times  load  levels  are  reached  or  exceeded  over  the  useful  life 
of  the  vehicle.  This  is  necessary  for  both  fatigue  design  criteria 
and  to  insure  an  acceptably  low  probability  of  gust  loads  in  excess 
of  structural  yield  limits.  This  load  exceedance  information  can 
be  obtained  by  the  use  of  power  spectral  techniques  for  the  repre- 
sentation of  atmospheric  turbulence  and  the  resultant  aircraft 
load  and  motion  responses . 

Several  different  models  have  been  developed  and  proposed  to 
represent  the  aircraft  gust  load  and  motion  responses.  Of  partic- 
ular interest  are  the  one-  and  two -degree -of- freedom  models  devel- 
oped by  Houbolt  (Refs.  1,2)  and  the  two-degree-of -freedom  models 
by  Peele  (Ref.  3).  Houbolt's  one-degree-of-freedom  (IDF)  models 
and  Peele 's  two-degree-of-freedom  (2DF)  models  were  used  to  esti- 
mate gust  loads  on  ten  aircraft  representative  of  those  subject  to 
Federal  Aviation  Regulation  23,  at  various  altitudes,  airspeeds, 
and  flight  configurations  in  Reference  A.  The  purpose  of  the 
study  reported  herein  is  to: 

1.  Apply  Houbolt's  longitudinal  and  lateral  2DF  models  to 
the  same  aircraft  and  flight  conditions. 

2.  Quantitatively  evaluate  the  three  techniques  in  order  to 
provide  data  for  possible  changes  in  the  Federal  Aviation 
Regulations  concerning  gust  design  procedures. 

3.  Develop  a handbook  implementing  a simplified  power 
spectral  gust  design  procedure. 


The  major  conclusions  of  this  study  are: 

1.  The  one-degree-of-freedom  Houbolt  method  is  recommended 
for  gust  load  estimates  if  the  method  must  be  simple 
enough  for  handbook- type  calculations. 

2.  Peele's  method  significantly  overestimates  the  gust 
loading  of  the  wing  on  the  aircraft  investigated. 

3.  The  two-degree-of-freedom  method  used  in  this  study  or 
one  of  its  variants  is  a useful  design  tool  for  the 
small  aircraft  manufacturer  with  access  to  a digital 
computer . 


2.  DESCRIPTION  OF  OUST  LOAD  ESTIMATION  TECHNIQUE 


2.1  Review  of  Load  Level  Exceedance  Estimation  Based  on 

Structural  Response  Parameter,  A , and  Rate  of  Zero 

Crossings,  Nq 

As  stated  in  the  introduction,  the  continuous  nature  of  turbu- 
lence produces  a distribution  of  aerodynamic  loads  in  time  which 
can  be  characterized  by  the  number  of  times  a load  level  is  ex- 
ceeded per  flight  hour  or  over  the  anticipated  life  of  the  air- 
craft. A suggested  use  of  this  load  exceedance  information 
(Ref.  1)  is  to  set  the  maximum  design  load  such  that  it  is  not 
exceeded  more  than  ten  times  in  30,000  hours  of  flight.  It  is  not 
within  the  scope  of  this  study  to  develop  or  evaluate  such  design 
load  criteria,  but  merely  to  point  out  that  this  approach  is  more 
realistic  than  the  discrete  gust  design  approach  required  in  FAR  23. 

The  primary  goal  of  this  study  is  the  comparative  evaluation 
of  three  different  methods  of  obtaining  these  load  exceedance  rates 
for  each  lifting  surface.  Therefore,  a brief  review  is  given  of 
the  general  manner  in  which  the  structural  response  parameter  A , 
and  rate  of  zerc  crossings,  No  , are  calculated,  their  physical 
significance,  ard  how  these  are  used  to  specify  the  load  exceedance 
rate.  This  is  then  followed  by  a description  of  Houbolt's  longi- 
tudinal and  lateral  2DF  gust  response  models. 

The  symbols  used  to  describe  aircraft  loads,  moments,  motions, 
and  their  sign  conventions  are  shown  in  Figure  2.1.  The  basic 
physics  of  the  i DF  models  used  are:  1)  that  the  aircraft  responds 
with  vertical  motion,  z , pitch  motion,  6 , and  horizontal  tail 
load,  P^  , to  a continuously  varying  vertical  component  of  atmo- 
spheric turbulent  velocity  w^  (Fig.  2.2),  and  2)  the  aircraft 
responds  with  side  motion,  y , yawing  motion,  ip  , and  lateral 
tail  load,  P^  , to  a continuously  varying  lateral  component  of 
turbulent  velocity,  v (Fig.  2.3).  The  basic  limitations  of 
this  approach  are  that  the  effect  of  the  turbulent  velocity  com- 
ponent in  the  forward  direction,  Ug  , is  neglected,  the  lateral 


y,  S,  q,  Y,  M 


*.  <£,P,X,L 


t,  r,  N 


Figure  2.1  Conventions  for  positive  position,  angular  rate, 
force  and  moment  about  each  axis  (moment  and 
angular  rate  follows  right  hand  rule) 

* The  position  and  force  are  defined  positive 
upwards  in  keeping  with  the  work  of  Ref.  2. 
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aircraft  response  does  not  include  roll,  and  the  longitudinal  and 
lateral  responses  are  uncoupled.  The  restrictions  are  justified 
on  the  basis  that  the  primary  aerodynamics  loading  of  the  lifting 


surfaces  and  fuselage  are  due  to  w and 

6 g 

roll  will  provide  conservative  estimates. 


and  that  neglecting 


An  additional  assumption  inherent  in  all  that  follows  is  that 


the  structural  response  modes  of  the  aircraft  are  well  damped, 
have  natural  frequencies  well  above  the  aerodynamic  response  modes, 
and  can  therefore  be  neglected  in  computing  the  aircraft  motion  and 
load  response  due  to  gusts.  This  is  usually  the  case  for  the  air- 
craft covered  by  FAR  23  due  to  their  relatively  small  size. 

Given  the  assumptions  listed  above,  it  is  possible  to  develop 
models  of  aircraft  responses  z » 0 ’ Pw  * PHT  to  w^  ; and  y , 

^ to  v in  either  the  time  or  frequency  domain  (where  ?w  , 

p 

HT  , PyT  are  the  loads  on  the  wing,  horizontal  tail,  vertical 
tail,  respectively).  Because  of  the  continuous  and  statistical 
nature  of  turbulence,  it  is  more  efficient  to  characterize  atmos- 
pheric turbulence  velocity  components  by  their  power  spectra 

$ (w)  , <p  (cj)  . This  then  requires  that  the  aircraft  gust 
g g 

response  functions  also  be  expressed  in  transfer  function  form 

Vu> 

H.(U> 

Hp  (u> 

HT 

Gust  load  estimates  are  then  obtained  from  each  of  the  output  power 
spectra  in  terms  of  the  input  power  spectra  and  the  transfer 
functions  shown  above 


‘ 


4>;-  (uj ) = H " (u) ) H-  (- w) <p  (w) 

Z Z Z g 


4>p  (<*>)  - Hp  (w)  Hp  ( - uj ) 4>  (u>) 

rHT  rHT  HT  g 


4>^.  Cto>  * H- (u>)  H- (-uj)  4>  (w) 

y y y o 


4>p  (u)  “ Hp  (u)  Hp  (>10)  (uj) 

vT  rVT  rVT  g 


(2.2) 


While  the  output  power  spectra  contain  all  of  the  necessary 

structural  loading  information,  it  is  not  in  a particularly  useful 

t hi 

form.  The  root  mean  square  value  of  the  i response  variable,  o^, 
is  of  more  interest  for  structural  design  purposes. 


■ (r)2 


(k) 


o..  1/2  rrQ  "l1'2 

(k)  dk 

8 J0  y J 


(2.3) 


These  integrals  can  be  written  in  terms  of  the  dimensional  frequency, 
uj  , or  its  nondimensional  equivalent,  k . 

k = fu  (2,A) 

where  c is  a reference  chord  length  i*nd  U is  forward  flight  speed. 
They  are  shown  here  in  terms  of  k , since  they  are  computed  in 
this  form  in  t ^ study. 

It  should  also  be  noted  that  the  upper  limit  on  these  inte- 
grals has  been  changed  from  infinity  to  a value  for  compu- 

tational purposes.  The  particular  value  of  cutoff  frequency 
chosen  and  its  effect  on  A and  Np  is  discussed  in  References  1 
and  2.  The  ratio  between  the  rms  values  of  the  output  and 

input  o is  defined  as  the  structural  response  parameter  A^n 
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substituting  (2.2)  into  (2.3)  yields 


An 


r-  k 

/' 


4>  (k> 


g 


-^1/2  „KC  ~ » -»l/2 


dk 


f ii 


001  V(k) 


Lro  g 


dk, 


w_ 


g 


w 


g 


AAn°w_ 


(2.5) 


Similarly  for  the  other  loading  variables 


HT 


r -kc  |Hp  (k)|2  ^ (k)  ^1/2 

^ 1 dk 


I 


g 


w 


w. 


g 


AHTc 


HT 


g 


a 


AX, 


°VT  = 


|H  (k) 
2’ 


2 <}>v  (k) 
g 


dk 


g 


1/2 


g 


A.  „0 

AX,  v 


These  structural  response  parameters,  A^ , where  i = An  , AX  , HT  , VT 
give  the  designer  a feel  for  the  magnitude  of  each  major  structural 
component  due  to  turbulence . 

Another  parameter,  important  for  structural  fatigue  estimates 
is  the  frequency  with  which  any  particular  load  level  is  exceeded. 
The  number  of  exceedances  of  a particular  acceleration  level  , An-^ 
per  unit  time  due  to  heavy  turbulence  is  given  by  the  relation 


where 


N. 


An 


is  the  rate  of  zero  acceleration  crossings  in  the 


positive  direction. 

N is  obtained  from  the  relation 
°An 


„ 1 °An  . U_ 

O.  Tir  O.  TfC 

An  An 


/0  k2  *4n<k>dk 


t 

I 


LO 


*An(k)dk 


1/2 


(2.7) 


It  can  be  seen  from  Eq.  (2.6)  that  the  gust  intensity, 


'w 


g 


structural  load  parameter,  A ^ , and  zero  crossing  parameter 

N completely  specify  the  number  of  times  any  desired  load 
o An 

factor  is  exceeded  per  unit  time.  This  same  relation  also  applies 


for  each  of  the  other  variables  of  interest  A Z 


and  P, 


HT  ’ VT 

Comparisons  of  the  gust  load  estimates  by  three  different  aircraft 
gust  response  models  are  therefore  presented  in  terns  of  the  pair 


of  parameters,  A 
listed  above. 


N 


associated  with  each  of  the  variables 


The  following  subsection  gives  the  characteristics  of  the 
atmospheric  turbulence  spectra  <J>__  (k)  , 6..  (k)  used  in  this 
analysis . 


4>  (k) 

Yw  ' 

g 


4>v  00 

g 


2.2  Model  of  Atmospheric  Turbulence 

For  all  three  load  estimation  techniques  examined  in  this 
study,  atmospheric  turbulence  is  characterized  by  the  von  Karman 
spectrum  for  both  vertical  and  lateral  gust  velocities 


2 1 + | (1.339U02 

<}>  (si)  = a — -= 

Twv  ' w_  ir  r 


g 


1 + (1.339L8) 


where 


> 


IT7^ 


(2.8) 


ft 


w 


= spatial  frequency  (rad/ft) 

= root  mean  squared  value  of  gust  velocity  (also  called 
g gust  intensity1) 

L = scale  of  turbulence 
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For  an  aircraft  traversing  a patch  of  turbulence,  the  spatial 
distribution  appears  as  a time  varying  function.  For  an  aircraft 
flying  at  velocity  U with  reference  wing  semichord  c/2  the 
gust  input  power  spectrum  can  be  expressed  in  terms  of  the  non- 
dimensional  frequency 


i fie  _ U)C 

k " T ‘ 2u 


(2.9) 


where 


to  = temporal  frequency  (rad/sec) 

, ...  22L1+T  U-339  §ik)2 

V(k)  = °w  — - 


S 


1 + (1.339  k)2 


TTT5 


(2.10) 


This  same  spectral  form  and  turbulence  scale  are  used  for  the 
lateral  gust  velocity,  v . In  addition,  the  gust  velocities  are 
assumed  to  be  uncorrelated.  The  basic  spectral  representation 
plus  the  assumption  concerning  relative  turbulence  scales  and 
lack  of  correlation  are  reasonably  well  justified  by  experimental 
data,  at  least  above  an  altitude  of  1000  ft.  The  scale  of  atmos- 
pheric turbulence  varies  with  altitude,  below  1000  ft,  atmospheric 
stability,  wind  shear  and  other  factors.  Estimates  of  turbulence 
scale  run  from  less  than  one  hundred  to  several  thousand  feet. 

For  the  purpose  of  this  study,  a value  of  L = 2500  ft  was  chosen 
for  most  analyses.  This  choice  is  not  crucial  since,  as  shown  in 
Section  3.4,  the  effect  of  scale  change  on  the  structural  loading 
parameter  A is  satisfactorily  given  by  the  relation 

a/3 


A(Ll) 


(2.11) 


for  scale  lengths  greater  than  750  ft. 


2.3  Model  of  Aircraft  Response  to  Gusts 

In  this  section  the  rationale  for  estimating  the  aerodynamic 
forces  on  the  lifting  surfaces  and  aircraft  fuselage  due  to  atmos- 
pheric turbulence  is  reviewed.  First,  the  basic  potential  flow 


j 


i 


model  is  described.  Then  the  influence  coefficient  relating  the 
aerodynamic  load  at  one  station  to  downwash  velocity  at  another 
station  is  developed.  Finally,  the  aerodynamic  loads  and  re- 
sulting moments  are  incorporated  in  the  two-degree-of-freedom 
longitudinal  and  lateral  equations  of  motion  in  order  to  calculate 
gust  response  transfer  functions,  structural  response  parameter, 

A , and  zero  crossing  parameters,  NQ  . In  addition,  equations 
for  the  equivalent  stability  derivatives  obtained  from  this  method 
are  developed  for  comparison  with  the  derivatives  used  in  Peele's 
methods  in  Reference  3.  These  equivalent  stability  derivatives 
should  also  be  useful  for  estimating  the  inaccuracy  of  assuming 
real  and  constant  values  for  these  functions  when  extracting  them 
from  flight  test  data  using  high  frequency  control  motion  input 
signals . 

The  basic  premise  of  Houbolt's  technique  is  that  the  loads 
and  velocity  flow  field  resulting  from  gust  velocity  can  be 
adequately  described  by  potential  flow  theory.  This  will  not  be 
true  if  the  flow  separates  from  a lifting  surface;  i.e.,  if  the 
aircraft  stalls  or  spins.  Nor  is  it  likely  to  accurately  account 
for  wing-body  interference  effects.  However,  for  the  prediction 
of  loads  on  the  major  surfaces  during  turbulence  encounters  at  or 
near  cruise  conditions  (low  angle  of  attack  and  sideslip)  it  has 
great  promise.  The  potential  flow  model  is  a series  of 
dipole  sources  which  produce  a line  load  as  illustrated  in 
Figure  2.4.  The  wing  load  is  approximated  by  two  line  loads,  the 
fuselage  by  four,  and  the  tail  surfaces  by  one  each.  This  figure 
also  shows  the  location  of  the  points  where  the  downwash  velocities 
are  computed  as  well  as  the  motion  variables  of  interest.  The  one- 
and  two-line  load  approximations  were  shown  in  Reference  2 to 
adequately  estimate  the  lift  curve  slope  of  an  elliptically  loaded 
wing  with  aspect  ratio  of  eight. 

The  basic  relations  obtained  from  the  line  load  model  is  the 
downwash  velocity  at  station  m due  to  a load  at  station  n . 

This  is  given  in  terms  of  the  influence  functions  , 3^ 
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(2.12) 


wm  iTpU  A ^amn  + mn ^ 


Where  for  the  downwash  station  ahead  of  the  load 

iks 


a + iB  = e 
mn  mn 


» •< 

Cl^so^  + iSl(so) 
— 


(2.13) 


and  for  the  downwash  station  downstream  of  the  load 

iks 


a + iB  = -e 
mn  mn 


where 


'°  [°o  + W + lSl<so>]  <2-14) 


s is  the  nondimens ional  distance  between  load  and  down- 
o 


wash  velocity 


so  = ‘ XDWNWSH  " xLOAdI2/,c 


(2.15) 


and  the  C and  S functions  are: 

/ 

cos  ks 


C = TTk  + 2 
o 


( 


1 - 


X/c 


/s2  + (X/c)2 


ds 


(2.16) 


oc 

:l<so>  -f 

so 

w -f 


X/c  cos  ks 
s2/s2  + (t/c) 


X/c  sin  ks 

s2/s2  + ( /c)2 


ds 


(2.17) 


ds 


(2.18) 


where 

s is  the  nondimensional  distance  variable,  positive  in 
the  downstream  direction 

s = - ~ (2.19) 

Both  k and  s are  a function  of  the  reference  chord,  c , 
appropriate  to  the  fuse] ape  or  particular  aerodynamic  surface  in 


question.  For  each  surface  the  reference  chord,  c , spanwise 
length  of  line  load,  A and  lift  curve  slope,  a , are  determined 
as  follows.  The  lift  curve  slope  is  a function  of  aspect  ratio, 

AR  , sweep  angle,  A , and  Mach  number,  M 


n AR 

1 + Vi  + (rsrn/O-  - ”2  cos2  A) 


(2.20) 


(2.21) 


where 

S 


c - | (2.22) 

is  the  area  of  the  aerodynamic  surface  in  question. 


A different  approach  is  used  for  the  calculation  of  A and  c 
for  the  fuselage. 

Af  = .707  kdf  (2.23) 

c f « .278  £ f (2.24) 

where 

subscript  f means  fuselage 
* length  of  fuselage 
d^  * average  diameter  of  fuselage 
K = constant  set  at  .75  to  produce  correct  subsonic 
aerodynamic  yawing  moment 

In  the  program  for  computing  gust  load  spectra,  the  wing, 
tail, and  fuselage  loads  are  all  calculated  together.  This  is  most 
easily  accomplished  with  all  equations  written  in  terms  of  the 
reduced  frequency  associated  with  the  wing,  . In  this  case 

the  a and  8 associated  for  example  with  the  horizontal  tail 
become 
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Similar  relations  also  apply  for  the  vertical  tail  and  the  fuselage. 

It  should  be  noted  that  these  equations  imply  that  not  only  is 
the  magnitude  of  the  downwash  velocity-load  relation  a function 
of  frequency  but  also  the  phase  angle  change  in  the  relationship 
of  these  variables . This  is  unique  to  this  modeling  technique  and 
should  therefore  provide  more  realism  than  Peele's  method.  In 
particular,  Peele's  method  relates  the  downwash  velocity  to  load  as 


% «=  PnKpe_ak  (2.29) 

where 

a is  a real  valued  lift  attenuation  factor  which  is  a 
function  of  Mach  number  and  sweep  angle  given  in 
Reference  3. 

Kp  constant  of  proportionality  determined  by  aircraft 
stability  derivatives 


X 1 

4-, 
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The  quantitative  significance  of  including  the  phase  angle 
variation  will  be  shown  in  Section  3.1. 

The  relations  presented  so  far  apply  to  both  the  longitudinal 
and  lateral  gust  load  responses.  At  this  point  the  development  of 
these  models  is  separated  and  the  longitudinal  response  model  is 
described  first,  followed  by  the  lateral  response  model. 


Referring  to  Figure  2.5  for  the  longitudinal  load-velocity 
control  points,  the  three  loads  produce  downwash  velocities  w^ 
at  the  control  points 

W1  ■ 7 C(“ll  + iBU)Pw1  + <a12  + l612)Pw23 
w2  ‘ y C(a21  + i621,Pw1  + + 1622)Pw23 


w. 


HT 


y C (a31  + iB31)Pw1  + (a32  + i032)Pw. 


+ ^ (a 


HT 


33  + i633)PHT^ 


(2.30) 


where 
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6 


22 
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The  downwash  at  the  wing  due  to  the  tail  load  is  neglected. 

The  downwash  velocities  must  also  satisfy  kinematic  contraints 

• • 

in  terms  of  the  vertical  aircraft  velocity,  z , pitch  rate,  6 

iwt 

and  gust  velocity  magnitude,  w . (The  oscillatory  term  e is 

O 

not  written  in  order  to  keep  the  notation  manageable,  but  is  under- 
stood to  be  part  of  all  input  and  output  variables.) 


c • 

w,  = z + (e  - 0^)0  - U0  - w 
1 V w 8 / g 
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(2.31) 
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The  aerodynamic  loads  can  be  expressed  in  terms  of  the  motion 
variables  z , 6 and  the  input  gust  velocity  w by  combining 


Eq.  (2 

.30) 

and  Eq . 

(2.31)  . 
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(2.32) 


The  longitudinal  aircraft  heaving  and  pitching  motions  are  des- 
cribed by  the  lift  and  pitching  moment  equations. 


mz  = P + P + p 

w^  v*2  HT 


mr 


Q 

;5  - \ (e» + if) + pw2  (ew  - h w)  - p 


HTeHT 


(2.33) 
(2  34) 


where  m is  aircraft  mass  and  r^  is  radius  of  gyration  about  the 
y axis. 

Equation  (2.32)  plus  the  two  equations  of  motion  form  a set 
from  which  the  motion  and  loads  can  be  determined.  This  set  of 
equations,  suitably  nondimensionalized , is  Houbolt's  longitudinal 
two-degree-of- freed' m gust  load  response  model  (Fig.  2.6). 

The  develoDmeot  of  the  lateral  load  gust  response  model  is 
quite  similar  to  the  longitudinal  case.  Figure  2.4b  illustrates 
the  locations  for  the  line  dipole  sour  :es  which  are  used  to  repre- 
sent the  forces  and  moments  produced  by  the  fuselage  and  vertical 
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Figure  2.6  Longitudinal  2DF  Gust  Load  Model 
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Kinematic  Equations  for  Sidewaah  Velocity  at  Control  Point* 
Vj  - y + (ef  + cf/2)i  - Ity  - Vg 

vj  - y + (e{  - Cf/2)i>  - U*  - 

v3  • y + (af  - 3/2cf)i>  - U*  * vg 

v4  * y + («f  - 5/2cf)i  - Ui*i  - vg 

y * (e„r  + cwr/2>^  * u*  " v« 


2.35) 


Sldewash  Velocity  - Side  Load  Equations 
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vvt  " ?7  sf  (avt  4 1Bvt> 


(2.36) 


Side  Force  Equation 


my  - 4 P^  + P^  + Pf^  + Pvt 


(2.37) 


y awing  Moment  Equation 

I„i  - <*f  ♦ cf)Pf  4 afPf  4 (ef  - cf)Pf3  4 (af  - 2cf)P£<(  4 .^P^ 

(2.38) 
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tail  due  to  lateral  gusts.  The  lateral  load  and  response  variables 
are  depicted  in  Fig.  2.7.  The  major  assumptions  in  this  model  are: 
1)  that  the  wing  contribution  to  lateral  forces  is  neglected;  2) 
aerodynamic  loads  on  the  body  do  not  significantly  contribute  to 
sidewash  on  the  vertical  tail  and  3)  aerodynamic  loads  on  the  ver- 
tical tail  do  not  significantly  contribute  to  the  sidewash  on  the 
fuselage.  Another  frequently  used  second-order  model  of  lateral 
motion  is  a combination  of  roll  and  yaw  with  sideforce  neglected. 
This  produces  the  familiar  Dutch  roll  motion  which  is  of  interest 
for  the  evaluation  of  lateral  handling  qualities.  This  model, 
however,  is  computationally  more  involved  since  wing  loads  must  be 
included.  In  addition,  it  is  not  as  close  to  Peele's  lateral  model 
so  that  comparisons  with  his  method  would  not  be  as  meaningful. 
Therefore,  while  the  Dutch  roll  model  is  of  interest  for  future 
studies,  it  will  not  be  used  in  this  study.  Proceeding  as  in  the 
longitudinal  case,  equations  relating  sidewash  velocity  to  aero- 
dynamic loads,  Eq.  (2.35),  and  to  kinematic  constraints,  Eq.  (2.36) 
are  written  along  with  the  side  force  and  moment  equations  (2.37) 
and  (2.38).  In  addition  to  these  relations,  slender  body  theory 
dictates  that  the  wake  potential  downstream  of  a body  of  revolution 
(the  fuselage)  should  be  zero.  In  order  to  accomplish  this,  the 
four  fuselage  loads  are  modified  as  explained  in  Appendix  A.  The 
side  force  and  moment  equations  become 
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This  set  of  seven  equations  suitably  nondimensionalized  (Fig. 
2.8)  forms  the  lateral  two -degree-of-fieedom  gust  load  response 
model  used  to  analyze  the  gusf  It  '1  exceedance  parameters  for  ten 
small  aircraft  as  discussed  in  the  next  section. 


Figure  2.8  Lateral  2DF  Gust  Load  Model 


This  equation  can  be  solved  numerically  for  the  line  loads  by 
inverting  A . 
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The  longitudinal  stability  derivatives  can  be  written  in  terms 
of  the  line  loads 
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The  line  load  derivatives  required  in  Eqs . (2.47)  - (2.50) 
are  proportional  to  the  coefficients  of  the  top  two  rows  of  matrix 
[C],  Eq.  (2.46). 
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Note  that  there  are  additional  pseudo  stability  derivatives 
in  this  model  due  to  the  change  in  the  magnitude  of  the  gust  velo- 
city along  the  x axis  of  the  vehicle. 
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For  turbulence  wavelengths  long  compared  to  the  length  of  the 
vehicle,  these  terms  should  be  small. 


The  equivalent  lateral  stability  derivates  C , CN  , 

ye  6 

C , and  CM  are  obtained  by  the  same  procedure.  The  yaw  rate,  r 
Nr  „ 

and  angle  of  side  slip,  6 , are  substituted  into  Eq . (2.35)  for  y 

and  v . This  new  equation  is  combined  with  Eq . (2.36),  solved 

numerically  for  the  line  load  derivatives  9Pf  / 3 B , etc.  , which 

are  combined  to  obtain  the  lateral  stability  derivatives. 

Numerical,  results  using  the.se  equations  aro  ?iven  in  Section 

3.5. 
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3.  GUST  RESPONSE  OF  TEN  GENERAL  AVIATION  AIRCRAFT 


The  mathematical  models  of  atmospheric  turbulence,  aircraft 
gust  response,  and  the  calculation  of  the  resulting  output  power 
spectra,  and  load  exceedance  parameters  A and  N0  were  coded 
in  Fortran  IV  for  solution  on  A. R. A. P.'s  PDP  11/70  digital  com- 
puter. The  program  was  run  for  ten  general  aviation  type  aircraft 
in  different  configurations  at  cruise  airspeed  and  several  altitudes 
(a  total  of  69  cases) . The  aircraft  and  conditions  chosen 
(Table  3.1)  include  all  of  the  cases  analyzed  in  Reference  4 by 
Peele's  2DF,  Houbolt's  IDF,  and  the  discrete  gust  methods.  The 
results  of  these  calculations  are  presented  primarily  in  terms  of 
pairs  of  load  exceedance  parameters 
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as  a function  of  altitude.  Additional  results  are  given  in  terms 
of  frequency  response  functions  which  are  the  square  of  the  ampli- 
tude of  the  transfer  functions  previously  discussed  (Section  2.1). 
The  frequency  response  functions  are  defined  similar  to  those  used 
in  References  2 and  5. 
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The  response  function  results  are  useful  for  determing  the  approxi- 
mate second-order  natural  frequency  and  damping  and  for  aid  in 
understanding  changes  in  the  load  exceedance  parameters.  These 
have,  for  the  most  part,  been  relegated  to  Appendix  B since  they  are 
auxilary  in  nature. 


The  load  exceedance  parameters  are  shown  as  function  of  alti- 
tude for  each  of  ten  aircraft  at  the  specified  flight  configurations 
(Figs.  3.1  - 3.10).  Two  basic  trends  are  seen  in  the  Houbolt  2DF 
results  contained  in  these  figures:  1)  a decrease  in  A^n  and, 
usually,  an  increase  in  Nq  with  altitude.  This  is  explained 
by  referring  to  a typical  se?  of  frequency  functions  as  shown  in 
Fig.  3.12a. 


For  airplane  #1  the  amplitude  squared  response  function,  fR  , 

increases  from  sea  level  to  20K  ft  by  about  30%  (Fig.  14a).  However, 

2 2 

this  must  be  multiplied  by  (irpUS/w)  to  convert  to  | An/w  | . The 

O 

density  terms  causes  the  load  factor  transfer  function  to  decrease 
with  altitude  thus  resulting  in  lower  values  of  A^  . Again, 
referring  to  Fig.  3.12a,  it  is  seen  that  an  increase  in  altitude 
amplifies  F^  near  and  above  the  resonant  frequency  and  attenuates 
F at  low  frequencies.  Therefore,  Nq  which  is  proportional  to 
the  radius  of  gyration  of  this  responseA¥ncreases  with  altitude. 


Another  trend  evident  in  the  results  is  that  moving  the  eg 
aft  (Fig.  ,'J.ll)  increases  the  structural  loading  parameters  A^n  , 
A,,  , and  decreases  the  zero  crossing  parame* ers , NQ  and  Nq 

This  is  due  to  lowering  the  short  period  frequency  .-’ndA^omparableAC 
lateral  frequency.  Thr.s  shift  increases  the  amouvt  of  low  frequency 
energy  in  the  output  spectrum  more  tian  the  decrease  in  high  fre- 
quency energy. 
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LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  #1 

□ HOUBOLT  2 0 F 
O PEELE  DERIV.  (NASA  CR  I975) 

A PEELE  OERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOL  T'S  AT  k • 0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  #2 

□ HOUBOLT  2DF 
O PEELE  DERIV.  (NASA  CR  I97S) 

A PEELE  DERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOLT’S  AT  k*0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  # 3 

□ HOUBOLT  2DF 
O PEELE  DERIV.  (NASA  CR-1975) 

A PEELE  OERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOLT'S  AT  k « 0) 
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LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  #4 

□ HOUBOlT  2 D F 
O PEELE  OERIV.  (NASA  CR-I97S) 

A PEELE  OERIV.  (AIRCRAFT  MFR) 

O PEELE  OERIV.  (AGREE  WITH 
HOUBOLT'S  AT  k « 0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  # 5 

□ HOUeOLT  2 D F 
O PEELE  DERIV.  (NASA  CR1975) 

A PEELE  DERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  *ITH 
HOUBOLT'S  AT  k = 0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  # 6 

□ H0U80LT  2 OF 
O PE ELE  DERIV.  (NASA  CRI975) 

A PEELE  DERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOLT'S  AT  k « 0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  #7 

□ houbolt  2D f 

O PEELS  DERIV.  (NASA  CR  I975) 

A PEELE  DERIV.  (AIRCRAFT  MFR) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOLT'S  AT  h « 0) 


LOAD  EXCEEDANCE  PARAMETERS  vs  ALTITUDE 

AIRCRAFT  #10 

□ NOUBOLT  2 OF 
O PEELE  DERIV.  (NASA  CR-1975) 

O PEELE  DERIV.  (AGREE  WITH 
HOUBOLT’S  AT  k«0) 


OAD  RESPONSE  FUNCT 


Figure  3.1 2d 


COMPARISON  OF  LOAD  EXCEEDANCE  PARAMETERS  USiNG 
HOUBOLT ' S IDF  AND  2DF  METHODS 

o HOUBOLT  IDF 
A HOUBOLT  2 DF 


AIRCRAFT  *\ 


O .02  .04  0 12  3 


Figure  3.13 


COMPARISON  OF  TAIL  LOAD  EXCEEDANCE  PARA 
METERS  BY  2 DF  AND  ODF  METHODS 

AIRCRAFT  #l 

□ HOUBOLT  2DF 

O OOF  METHOO 

O PEELE  2 D F DERIV. 

NASA  CR-1975 
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COMPARISON  OF  TAIL  LOAD  EXCEEDANCE  PARA 
METERS  BY  2 DF  AND  ODF  METHODS 

AIRCRAFT  *5 

□ HOUBOLT  2DF 

O OOF  METHOD 

O PEELE  2 0 F DERIV. 

NASA  C R - 19  75 
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3.1  Comparison  of  Houbolt's  and  Peele's  2DF  Methods. 


Peele's  method  makes  use  of  the  standard  longitudinal  short 

period  model  and  lateral  snaking  model  with  stability  derivatives 

modified  by  exponential  lift  attenuation  factors.  The  results  of 
y 

this  approach  are  obviously  dependant  on  the  values  of  the  stability 
derivatives  chosen  to  represent  a particular  aircraft.  The  FAA 
study,  Ref.  4,  computes  the  load  exceedance  parameters  based  on 
derivatives;  1)  supplied  by  the  manufacturers,  and  2)  estimated 
using  methods  described  in  NASA  CR  1973  (Ref.  7). 

These  results,  also  included  on  Figs.  3.1  - 3.10,  illustrate 
a major  disadvantage  of  Peele's  method,  namely  the  significant  vari- 
ation in  the  results  due  to  variations  in  the  estimates  of  the  neces- 
sary stability  derivatives.  The  requirement  to  standardize  deriva- 
tive estimation  proceedures  becomes  evident  if  such  a technique  is 
to  be  incorporated  in  FAR  23.  This  in  itself  would  be  quite  an 
undertaking.  (It  should  be  noted  that  some  of  the  differences  in 

N estimates  taken  from  Ref.  4 are  due  to  an  inconsistent  choice 
o 

of  cut-off  frequency.) 

The  results  of  Houbolt's  and  Peele's  techniques  differ  both 
because  they  can  start  with  different  derivatives  at  zero  frequency 
and  because  the  functional  dependence  of  the  derivatives  on  fre- 
quency is  different.  To  separate  these  two  effects,  additional 
calculations  were  run  using  Peele'  method  with  the  derivatives 
set  equal  to  the  equivalent  zero  frequency  values  used  in  Houbolt's 
2DF  model.  These  results  are  also  shown  in  Figs.  3.1  - 3.10. 

The  basic  result  of  this  latter  comparison  is  that  in  all 

cases  values  of  A^n  and  N are  larger  when  computed  by 

°An 

Peele's  method.  This  is  due  to  Peele's  method  predicting  a larger 
normal  load  response  at  all  frequencies,  particularly  above  1 hz 
as  illustrated  in  Fig.  3.12. 

Comparisons  of  the  lateral  responses  indicate  that  values  of 
A.  are  smaller  and  N are  larger  when  computed  bv  Peele's 

ii  c O . , 

A C 

method.  This  also  is  explained  bv  referring  to  typical  transfer 
function  plots  (Fig.  3.12)  which  indicate  that  Peele's  method  under- 
estimates the  lateral  response  at  low  frequency  and  overeat imates 
in  the  high  frequency  range. 
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Figure  3.10a 


The  basic  conclusion  arising  from  these  comparisons  is  that , 
by  neglecting  the  phase  angle  shift  in  lift  attenuation.  Peele's 
method  overestimates  the  main  load  exceedance  parameter  A^n  . For 
this  reason  and  because  of  the  sensitivity  to  stability  derivative 
estimates  it  is  not  recommended  that  Peele's  method  be  incorporated 
into  FAR  23 . 

3.2  Comparison  of  Results  Obtained  Using  Houbolt's  IDF  and  2DF  Models 

Reference  1 describes  a IDF  model  with  which  the  load  exceedance 
parameters  can  be  calculated  manually  using  two  charts . In  tnis  pro- 

ceedure  the  aircraft  is  assumed  to  respond  longitudinally  in  heave 
(only  the  lift  equation  is  used  so  the  aircraft  is  characterized 

completely  by  the  density  parameter  p = 2w/apcgS  . This  method 
has  great  appeal  since  it  correctly  models  the  unsteady  lift  on  the 
wing  and  is  quite  simple  to  apply.  The  drawback  of  this  approach 
is  that  by  neglecting  the  pitch  response,  the  load  alleviating 
tendency  of  the  aircraft  to  weathercock  (hold  angle  of  attack  con- 
stant) is  not  included.  How  important  this  effect  is  for  FAR  23 
aircraft  is  illustrated  in  Fig.  3.13.  Here  the  normal  load  exceed- 
ance parameters,  calculated  using  Houbolt's  IDF  (Ref.  1)  and  2DF 
(Section  2)  methods,  are  presented.  For  the  sake  of  brevity  the 
comparison  is  limited  to  three  aircraft  chosen  to  cover  the  range 
of  FAR  23  aircraft.  They  are  a light  single  engine  aircraft,  a 
twenty-passenger  turboprop,  and  a small  jet  business  transport. 

This  comparison  of  results  indicates  that  the  inclusion  of 
the  pitch  response  reduces  the  nus  normal  load  factor  estimate  by 
as  much  as  507,,  ac  sea  level,  but  can  even  underestimate  at  high 
altitude.  This  is  because  the  damping  of  tne  short  period  mode 
decreases  with  altitude,  thus  increasing  the  resonant  response  of 
the  second-order  2DF  model.  This  is  a significant  difference  but 
it  should  be  noted  that  manual  calculation  of  the  load  exceedance 
parameter  by  Houbolt  's  2DF  method  using  cha~ts  is  not  feasible  duo 
to  the  large  number  of  necessary  input  parameter-  It  is,  therefore, 
recommended  thaL  the  IDF  method  be  used  if  manual  calculations  arc 
to  be  employed. 
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3.3  Comparison  of  Load  Exceedance  Parameters  Based  on  2DF  and  ODF 
Models . 

It  is  desirable  to  obtain  a simple  method  for  estimating  gust 
load  exceedance  parameters  for  the  vertical  and  horizontal  tail- 
surfaces.  The  loads  on  the  horizontal  tail  are  a strong  function 
of  the  unsteady  wing  downwash,  therefore,  they  can  not  realistically 
be  calculated  independent  of  wing  conditions.  The  vertical  tail, 
however,  is  somewhat  above  the  fuselage  and  less  strongly  influenced 
by  the  relatively  weak  fuselage  sidewash.  In  addition  to  weak  coup- 
ling with  the  flow,  induced  by  other  parte  of  the  aircraft,  the  gust 
load  on  the  vertical  tail  produces  sideways  accelerations  which  are  small 
relative  to  the  wing  induced  vertical  accelerations.  This  would 
indicate  that  the  attenuating  effect  of  the  sideways  motion  can  pos- 
sibly be  neglected.  If  the  aircraft  also  does  not  develop  signifi- 
cant yaw  accelerations  it  is  worth  considering  a ODF  model  for 
vertical  tail  gust  loads.  For  trial  purposes  this  is  equivalent 
to  a 2DF  model  with  infinite  mass  and  moment  of  inertia.  The 
resulting  load  exceedance  parameters , A and  Nq  , for  both  the 
vertical  and  horizontal  tail  loads  for  three  typical  FAR  23  aircraft 
are  shown  in  Fig.  3.14  and  compared  with  the  results  from  the  stan- 
dard 2DF  model. 

The  ODF  model  consistantly  overestimates  the  tail  loads  due 
to  the  overestimation  of  the  low  frequency  tail  load  response  as  shown 
in  Fig.  B. 6 of  appendix  B.  This  suggests  that  a IDF  response  model 
allowing  yaw  motion  would  probably  do  a better  job.  In  lieu  of  devel- 
oping such  a model,  the  ODF  model  is  recommended  for  the  manual  esti- 
mation of  gust  loads  on  vertical  tails. 

Estimates  of  A and  N from  Ret.  4 given  by  Peele’s  method 

o 

at  sea  level  are  also  shown.  Structural  loading  is  high  compared 
with  Houbolt’s  2DF  method  and  not  consistently  related  to  Houbolt's 
IDF  method. 

3.4  Effects  of  the  Scale  of  Atmospheric  Turbulence  on  Gust  Load 

Excppdance  Parameters 

The  actual  scale  of  atmospheric  turbulence  goes  iron-  less  than 
100  feet  (close  to  the  ground)  up  to  five  thousand  feet.  The  prac- 
tical range  of  interest  for  aircraft  design  purposes  is  roughly 
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750  ft  s L £ 2500  ft.  The  particular  value  chosen  does  have  a sig- 
nificant effect  on  the  loading  parameter  A . Fortunately,  in  the 
range  of  frequencies  contributing  most  to  the  output  power  spectra, 
the  von  Karman  input  power  spectrum  is  approximately 

~ -522  °a)2  ([?)  w 5/3  (3*7) 


The  root  mean  square  value  for  a response  variable,  such  as  An  , 
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or  approximately 
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a^n  and,  therefore,  the  structual  loading  parameter  A^n  vary 
inversely  with  the  cube  root  of  the  turbulence  scale,  L . 


AA„<L2> 


AA„<L1> 


(3.10) 


This  relationship  is  accurate  to  within  5%  for  all  of  the  FAR  23 
aircraft  and  flight  conditions  studied.  This  is  illustrated  in  Fig 
3.15  for  aircraft  number  six  which  with  its  high  cruise  velocity 
should  have  the  largest  errors. 

The  exact  and  approximate  values  of  A are  shown  using 
Houbolt's  2DF  and  IDF  models.  The  error  is  greater  with  the  IDF 
model  since  this  model  responds  more  to  low  frequency  turbulence 
than  the  2DF  model  does . 

The  number  of  aero  crossings , Nq  , is  proportional  to  the 

An 

ratio  of  and  o^n  . Since  both  of  these  quantities  are  approxi- 

mately proportional  to  (L)  ^73their  ratio  is  independent  of  scale. 
Based  on  the  calculations  made  in  this  study,  for  aircraft  numbers 
1,  5,  6,  and  10,  this  approximation  is  957»  accurate  for  the  scale 
range  750  ft  < L < 2500  ft. 
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VARIATION  OF  STRUCTURAL  LOADING  PARAMETER 
Aa  , WITH  SCALE  OF  TURBULENCE,  L 


AIRCRAFT  w 6 
H = 40  KFT 
IAS  = 878  FPS 


O H0U30LT  IDF  MODEL 


2500 


2000 


The  two  conclusions  to  be  drawn  from  this  investigation  into 
scale  effects  for  FAR  23  aircraft  are: 

1)  Since  A varies  significantly  with  turbulence  scale,  a 
particular  value  of  L should  be  specified  for  design  purposes  or 
comparative  analyses . 

2)  The  particular  design  value  chosen  for  L is  not  crucial 

since  the  variations  of  A and  N with  scale  are  well  behaved  and 

o 

known . 

3.5  Stability  Derivative  Evaluation 

The  line  load  model  can  be  expressed  in  stability  derivative 
form  as  described  in  Section  2.3.  This  is  useful  since  the  equi- 
valent model  stability  derivatives  can  be  compared  with  those  ob- 
tained from  experimental  data,  Houbolt's  or  other  theoretical 
estimation  techniques,  and  the  accuracy  of  the  model  at  zero  fre- 
quency evaluated.  In  addition,  the  phase  angle  change  between  down- 
wash  velocity  and  induced  loads  gives  rise  to  an  imaginary  component 
in  each  derivative.  The  importance  of  this  phase  angle  effect  can 
be  estimated  by  the  relative  size  of  the  imaginary  and  real  components 
of  each  derivative. 

Figure  3.16  presents  the  complex  stability  derivatives  for  air- 
craft 1,  5,  and  6 at  cruise  conditions.  The  derivatives  are  eval- 
uated at  five  frequencies;  kw;  .001,  .05,  .1,  .25,  and  .5.  Values 
shown  at  the  lowest  frequencies  can  be  compared  with  the  steady  state 
stability  derivative  estimates  obtained  using  the  methods  given  in 
Refs.  7-12,  with  results  contained  in  Ref.  4.  Figure  3.16  shows 
derivative  values  obtained  using  NASA  CR  1975  (Ref.  7).  The  dif- 
ference between  derivative  values  obtained  from  Houbolt's  model  and 
the  technique  of  Ref.  7 is  the  same  magnitude  as  the  spread  in 
estimates  using  the  techniques  of  Refs.  7 - 12,  i.e.,  its  accuracy 
is  comparable  to  that  of  other  methods.  Also,  as  would  be  expected, 
differences  in  the  rotary  derivative  estimates  ( C,  , Cw  , C , 

q ‘q 

C*.  ) are  generally  larger  than  the  differences  in  the  static  de- 

r 

rivative  estimates. 
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Figure  3.16  also  indicates  that  all  of  the  derivatives  have 
significant  imaginary  components  when  the  aircraft  is  responding 
to  gust  or  control  inputs  at  frequencies  above  “0.5  In 
fact,  the  imaginary  components  of  Cmq  and  are  quite 

significant  for  > 0.05  . For  aircraft  1,  5,  and  6 this 
corresponds  to  dimensional  frequencies  of  .5,  .8,  and  1.25  Hz, 
respectively.  This  indicates  that  the  complex  representation  of 
these  derivatives  is  necessary  for  both  gust  load  modeling  and 
the  analysis  of  flight  test  data  with  gust  or  control  inputs  of 
frequencies  above  1 Hz. 


4.  SUMMARY  AND  CONCLUSIONS 


At  present  the  gust  load  estimates  specified  in  FAR  23  are 
based  on  encountering  a discrete  gust.  However,  commercial  and 
military  aircraft  use  more  sophisticated  power  spectral  design 
techniques  which  more  properly  account  for  the  continuous  and 
random  nature  of  turbulence.  Should  FAR  23  be  updated  to  include 
power  spectral  design  requirements  and,  if  so,  which  particular 
techniques  are  most  appropriate?  Answers  to  these  questions  must 
be  based  on  a knowledge  of: 

1.  The  accuracy  of  the  various  methods  in  predicting 
gust  load  power  spectra. 

2.  The  complexity  of  the  different  methods. 

3.  Ultimately,  how  well  does  each  method  discriminate, 
from  the  gust  load  design  viewpoint,  safe  aircraft 
designs  from  unsafe  ones. 

The  purpose  of  this  study  is  to  evaluate  the  relative  accu- 
racy and  complexity  of  three  power  spectral  techniques  and  to 
make  recommendations  based  on  these  evaluations.  It  is  believed 
that  these  recommendations  are  limited  by  the  fact  that  no  com- 
parison with  experimental  data  can  be  made  at  this  time. 

In  this  study,  response  functions  and  gust  load  exceedance 
parameters  have  been  calculated  for  ten  FAR  23  type  aircraft  at 
different  flight  conditions  and  configurations  using  Houbolt's 
2DF  technique.  These  results  have  been  compared  with  similar 
data  obtained  by  using  Peele's  2DF  method  and  Houbolt's  IDF 
method.  The  following  conclusions  are  given  based  on  these 
comparisons . 

Peele's  2DF  method,  while  somewhat  simpler  than  Houbolt's 
2DF  method,  requires  computerized  solutions  if  the  cutoff  fre- 
quency is  related  to  the  aspect  ratio  of  the  wing.  In  addition, 
this  method  neglects  the  significant  phase  angle  shift  with 
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frequency  between  aerodynamic  load  and  gust  velocity.  Since  the 
method  is  not  necessarily  more  accurate  than  Houbolt's  IDF 
procedure  and  it  is  certainly  move  complicated,  it  is  not  recom- 
mended for  use  in  gust  load  design  procedures. 

Houbolt's  IDF  method  overestimates  both  the  rms  gust  loading 
and  the  number  of  zero  crossings.  It  is,  however,  simple  to 
apply  and  is  the  best  candidate  for  updating  from  discrete  gust 
load  to  power  spectral  techniques. 

Houbolt's  2DF  method  is  based  on  geometric  and  inertial 
aircraft  characteristics.  These  are  mush  easier  inputs  to 
specify  than  the  stability  derivatives  or  frequency  and  damping 
characteristics  required  in  Peele’s  method.  On  the  other  hand, 
this  method  does  not  easily  allow  incorporation  of  flight  test 
data  results  into  the  model.  Because  of  the  large  number  of 
inputs  required,  manual  solution  with  the  aid  of  charts  is  not 
feasible.  Depending  on  the  approximations  made  in  determining 
the  influence  functions,  the  computer  code  can  fit  in  a program- 
mable calculation  (400  bytes)  or  require  21,000  bytes  of  storage. 
Further  refinements  in  the  2DF  models  follovjed  by  detailed  pro- 
gram documentation  are  necessary  before  this  method  can  become  a 
useful  design  tool. 

Houbolt's  2DF  model  appears  to  include  most  of  the  gust 
velocity-loading  phenomena  thought  to  be  significant,  although 
it  still  remains  to  be  validated  with  experimental  data. 


APPENDIX  A 


REMOVING  THE  FUSELAGE  WAKE 
WITH  A FIFTH  LOAD 

The  strength  of  the  wake  potential  due  to  a uniform  line 
load  of  intensity  p downsteram  a distance  x is 

ili)X 


= . JL 

pu 


For  n line  loads  a distance  Cf  apart  from  each  other  (see 
Fig.  2.7),  the  total  strength  of  the  wake  potential  is 


1 n 

1 P • e 
PU 


[x  + cf  (n- j ) ] 


where  x is  the  distance  between  any  desired  location  downstream 
of  the  aircraft  and  the  rearmost  line  load  pn  . For  a body  of 
revolution,  slender  body  theory  requires  that  the  strength  of  the 
wake  be  zero  downstream  of  the  body. 

A fictitious  5th  load  is  added  to  the  existing  four-panel 
representation  of  the  fuselage  to  satisfy  this  condition.  At 
the  point  of  this  added  load 
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The  5th  line  load  intensity  must,  therefore,  be 
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P5  = - l PJ' 
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Since  all  of  the  line  load  intensities  extend  over  the  same 
width  A , the  total  loads  p^  are  related  in  the  same  fashion 
as  the  intensities  p^  . 

The  additional  load,  p^  , is  included  in  the  side  force 
and  yawing  moment  equations  implicitly  in  terms  of  the  original 
four  loads  and  a moment  arm,  e^  - 3c^  . 
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APPENDIX  B 

This  appendix  contains  all  of  the  frequency  response  func- 
tions calculated  using  Houbolt's  2DF  method.  The  functions  are 

f 

useful  for  understanding  the  trends  seen  in  the  load  exceedance 
parameters  and  for  determining  each  aircraft's  approximate  lat- 
eral and  longitudinal  damping  and  natural  frequency.  The  response 
functions  relating  the  input  and  output  power  spectra  are  pro- 
portional to  the  square  of  the  amplitude  of  the  transfer  functions 
j 1 Eq.  (3.1)  and  Eq.  (3.6) . 
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